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Abstract 

The fragmentation functions for gluons to split into P-wave heavy quarkonium 
states are calculated to leading order in the QCD coupling constant. Long-distance 
effects are factored into two nonperturbative parameters: the derivative of the radial 
wavefunction at the origin and a second parameter related to the probability for a 
heavy-quark-antiquark pair that is produced in a color-octet S-wave state to form a 
color-singlet P-wave bound state. The fragmentation probabilities for a high transverse 
momentum gluon to split into the P-wave charmonium states XcO> Xd, and Xc2 are 
estimated to be 0.4 x 10 -4 , 1.8 x 10 -4 , and 2.4 x 10 -4 , respectively. This fragmentation 
process may account for a significant fraction of the rate for the inclusive production 
of Xcj a t large transverse momentum in pp colliders. 
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Heavy quarkonium plays an important role in high energy collider physics, because 
these states can probe physical processes at short distances of order I/itlq, where tuq is 
the heavy quark mass. Of particular importance experimentally are the 1 S-wave states 
of charmonium and bottomonium, which have very clean signatures through their leptonic 
decay modes, and the J ++ P-wave states with J = 0, 1, 2, which can also be observed 
through their radiative transitions into the 1 states. In most previous studies of the direct 
production of heavy quarkonium JlJ , the dominant production mechanisms were assumed to 
be given by the Feynman diagrams that were lowest order in the QCD coupling constant 
a s . It has recently been pointed out that the dominant mechanism for heavy quarkonium 
production at large transverse momentum px is fragmentation, the production of a high 
energy parton with even larger transverse momentum which subsequently decays into the 
quarkonium state plus other partons ||. While this mechanism is often of higher order in the 
QCD coupling constant a s than conventional mechanisms, it is enhanced at large transverse 
momentum by powers of PT/ m Q, and thus dominates at sufficiently large pt- 

The fragmentation of a parton % into any hadron H is described by a universal frag- 
mentation function Di_>jj(z, //), where z is the longitudinal momentum fraction of the hadron 
relative to the parton and fi is a factorization scale of order Pt • If the fragmentation 
function is known at some initial momentum scale [Iq, then it can be determined at larger 
momentum scales \x by solving the Altarelli-Parisi evolution equations which sum up the 
leading logarithms of [if [Iq. In Ref. 0, it was shown that in the case of heavy quarkonium, 
the fragmentation function D(z,rriQ) at an initial scale of order ijiq can be calculated using 
perturbative QCD. The initial fragmentation functions for gluons to split into S-wave states 
of heavy quarkonium were calculated to leading order in a s [0. The fragmentation func- 
tions for heavy quarks to split into S-wave states have also been calculated to leading order 
[|, §, |J, and these calculations have recently been extended to the P-wave states |]]. 

In this paper, we calculate the fragmentation functions for gluons to split into the 
P-wave states to leading order in a s . For the sake of clarity, we describe the calculation in 
terms of the lowest P-wave states of the charmonium system: the J PC = J ++ states Xcj, 
with J = 0, 1, 2, and the 1 H state h c . Our results apply equally well to the higher P-wave 
states of charmonium, as well as to the P-wave states of the bottomonium system. While a 
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gluon can split into Xcj at order a 2 through the Feynman diagram in Figure 1, gluon splitting 
into the h c occurs first at one order higher in a s . Since 1 H states of heavy quarkonium like 
the h c are difficult to observe experimentally, we concentrate in this paper on the J ++ states. 

In charmonium, the charmed quark and antiquark are nonrelativistic with typical 
velocity «<1 and typical separation l/(m c v). Our calculation of the fragmentation func- 
tion is based on separating short-distance effects involving the scale 1 / m c from long-distance 
effects involving scales of order l/(m c v) or larger. There are two distinct mechanisms that 
contribute to the fragmentation function at leading order in v ||, and we will refer to them 
as the color-singlet mechanism and the color- octet mechanism. The color-singlet mechanism 
is the production of a cc pair in a color-singlet 3 Pj state with separation of order l/m c in 
the quarkonium rest frame. The subsequent formation of the Xcj is a long-distance process 
with probability of order v 5 . In addition to the volume factor (m c v/m c ) 3 , there is an extra 
suppression factor of v 2 from the wavefunction of the P-state near the origin. The color- 
octet mechanism is the production of a cc pair in a color-octet 3 Si state with separation of 
order l/m c . The subsequent formation of the Xcj can proceed either through the dominant 
|cc) component of the Xcj wavefunction or through the small \ccg) component, which has a 
probability of order v 2 . In the first case, the cc pair must radiate a soft gluon to make a 
transition to the color-singlet 3 Pj |cc) state. In the second case, a soft gluon must combine 
with the cc pair to form a color-singlet \ccg) state. In either case, the probability is of order 
v 5 , with a volume factor of v 3 and an additional suppression of v 2 coming either from the 
probability of radiating a soft gluon or from the small probability of the \ccg) component 
of the wavefunction. Since the color-singlet mechanism and the color-octet mechanism con- 
tribute to the fragmentation function at the same order in v, they must both be included 
for a consistent calculation. 

Separating effects due to short distances of order 1 / m c from those of longer distance 
scales requires the introduction of an arbitrary factorization scale A in the range m c v <C 
A <C m c . The fragmentation functions for heavy quarkonium satisfy factorization formulas 
that involve this arbitrary scale. At leading order in v 2 , the factorization formula for the 
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fragmentation function for g — > Xcj has two terms: 

D^ XcJ (z,m c ) = ^ld[ J \z,A) + (2J + l)^-d 8 (z) , (1) 
m c m c 

where Hi and H%(A) are nonperturbative long-distance factors associated with the color- 
singlet and color-octet mechanisms, respectively. The short-distance factors d[ J \z,A) and 
ds(z) can be calculated using perturbation expansions in a s (m c ). They are proportional 
to the fragmentation functions for a gluon to split into a cc pair with vanishing relative 

momentum and in the appropriate color-spin-orbital state: color-singlet 3 Pj state for d[ J ^ and 
color-octet 3 Si state for d$. Note that in the factorization formula (|I|), the only dependence 

on A is in d[ J ^ and H' % . This simple form holds if the coefficients are calculated at most to 
next-to- leading order in a s . Beyond that order, d$ also acquires a weak dependence on A 
and J. 

The nonperturbative parameters Hi and ifg can be rigorously defined as matrix 
elements of 4-quark operators in nonrelativistic QCD ||. Their dependence on A is given by 
renormalization group equations whose coefficients can be calculated as perturbation series 
in a s (A) ||. To order a s , Hi is scale invariant and H' s satisfies 

If the factorization scale A is chosen to be much less than m c , this equation can be used to 
sum up large logarithms of m c /A: 

« = H'(A) + -HUg (^-] Hi , (3) 

27 (J \a s {m c ) ) 

where (3q = (33 — 2n/)/6 is the first coefficient in the beta function for QCD with rif flavors 
of light quarks. The parameter Hi can be related to the derivative of the nonrelativistic 
radial wavefunction at the origin for the P-wave states: 
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This parameter can be determined phenomenologically from the annihilation rates of the 
Xcj states. Using recent high precision measurements of the light hadronic decay rates of 
Xd and Xc2, Hi has been determined to be approximately 15 MeV [It]]. The parameter 



was introduced in Ref. |IT] in a calculation of the rate for the decay b — > Xcj + AT, which 
also receives contributions from both the color-singlet and color-octet mechanisms for Xcj 
production. The prime on H' % is a reminder that this parameter is not related in any rigorous 
way to the corresponding parameter H 8 that appears in decays of the Xcj states into light 
hadrons. Using data on the inclusive decays of S-mesons into charmonium, its value for 
A = m c has been estimated to be ifg(m c ) « 3 MeV JIT). This parameter also enters into the 



inclusive decay rate of the T into P-wave charmonium states |12 



We now turn to the calculation of the coefficient d[ J \z, A) in the color-singlet contri- 
bution to the fragmentation function. We follow the method and notation of Ref. [Q]. Let 
A a denote the amplitude for g* — > cc( 3 Pi, 1) + g corresponding to the Feynman diagram 
in Figure 1. The cc pair have almost equal momenta, and are in a color-singlet 3 Pj state. 
The amplitude A a can be written down in terms of R' P {0) using standard Feynman rules 
for quarkonium processes O]. Multiplying A a by its complex conjugate and summing over 
final colors and spins, we obtain the following generic form 

Y^,AxA*p = — ( Aj(s)(-g a/3 ) + Bj(s)p a pp + C,j{s){p a qp + q a pp) + Dj(s)q a qf3 ) , (5) 
trie \ J 

where p and q are the 4- momenta of the cc pair and the fragmenting gluon g*, respectively, 
and s = q 2 . The strategy is to reduce this expression in the limit go 3> m c to the polarization 
sum (—g a /3 + •••) for an on-shell gluon multiplied by a function of s and z, where z is the 
longitudinal momentum fraction of the cc pair relative to the fragmenting gluon. Terms in 
(ID that are proportional to q a or q^ can be dropped, because in axial gauge, q a and qp are 
of order m^/g when contracted with the numerator of the propagator of the virtual gluon. 
In the PaPp term, we can set p = zq + p±, where p± is the transverse part of the 4-vector p. 
After averaging over the directions of p±, p a pp can be replaced by —g a ppl/2, up to terms 
that are suppressed in axial gauge. The terms in (|5|) that contribute to fragmentation then 
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reduce to 



-gap) 



(6) 



Energy-momentum conservation in the form s = (p]_ + Am*)/ z +p±/(l — z) can be used to 
eliminate p]_ in (||) in favor of s and 2. The fragmentation probability is obtained by dividing 
the coefficient of (—g a p) by s 2 for the propagator of the virtual gluon, and then integrating 
over the phase space of the cc pair and the gluon in the final state. The phase space integral 
can be expressed compactly in terms of integrals over s and z 0]. The resulting expression 

for the integral over z of d[ J \z,A) is 



dzd[ J) (z,A) = 

16tH A min (A) J 4m*/ 



1 r°° 1 r 1 1 1 (a ,\ (! -4m 2 ) _ , A ff7 . 

ds I dz — \ Aj(s) + ± ^- -£j(s) . (7) 



We have anticipated the presence of an infrared divergence associated with a soft gluon in 
the final state by imposing a lower cutoff A on the energy of the gluon in the quarkonium rest 
frame. This translates into a lower limit on s: s m j n (A) = 4m 2 (l + A/m c ). The calculations 
of the functions Aj(s) and Bj(s) in (0) involves some rather complicated algebra, but the 
final results are relatively simple. Interchanging orders of integration in (|7|), we can read off 

the functions d[ J \z, A): 



d[ J \z,A) 



ds 



27 J Ami, '* s 2 (s-4m 2 c y 



fj(s,z) , z < 1 + 



ds 



nr. 



27J Sinin (A) s 2 (s-4m 2 ) 4 



A 

m c 



fj(s,z) , 2 > ( 1 + — 

m c 



(8) 
(9) 



where 



/o(s,*) 
fi(s,z) 

h{s,z) 



[s — 12m 



2\2 



s-4m 2 ) 2 - 2(1 - z)(zs - 4m 2 )s 



6 s 
2 



(s-4m 2 ) 2 - 2(1 - z)(zs -4m 2 )(s 



> - 4m 2 ) 2 (s 2 + 96m*) - 2(1 - z)(zs - 4m 2 )s(s 2 - 24sm 2 + 96m 4 c ] 



(10) 

(11) 
(.12) 



For z < (1 + A/m c ) 1 , the integral over s in (|8|) can be calculated straightforwardly. The 
cutoff A can be set to zero everywhere except in terms proportional to 1/(1 — z), which 



diverge upon integrating over z. In the 1/(1 — 2) terms, the limit A <C m c must be taken 
more carefully, and it gives rise to a plus-distribution: 

' e(l-z-— ) — ► — log— 5(1 -z). (13) 



1 — z V m c / (1 — z)+ m c 

For z > (1 + A/m c ) _1 , the limit A <C m c can be taken only after evaluating the integral over 
s in @. This gives rise to additional endpoint contributions proportional to 8(1 — z). Our 
final result for the short distance factor multiplying Hi/m c in the fragmentation function is 



4 j, (,,a) = 



(2.7+1) - - ■ + (Qj-(2J+l)log—)s(l-z) + Pj(z) 
1 1 — z)+ \ m c / 



(14) 



where the coefficients Qj are 
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Qo = g , Qi = f , Q 2 = ^ , (15) 



and the functions Pj(z) are 



p oW = £2i_^£) ^^j, (16) 

p i(z) . _M±M, (17) 

= + 9(2 - z) log(l - z) . (18) 



We next consider the color-octet coefficient d%(z) in the fragmentation formula (p]). 
At leading order in a s , this contribution to the fragmentation function comes from the 
subprocess g* —>■ cc( 3 S*i,8) given by the Feynman diagram in Figure 2. The c and c have 
equal momenta q/2, and are in a color-octet 3 Si state. The projection onto this state can 
be reduced to a simple Feynman rule: 

v(q/2)u(q/2) -> -^Ll^ ?(q)(fi+2m c ) , (19) 
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where e fl (q) is the polarization 4- vector of the 3 Si state and i, j, and a are the color indices of 
the quark, antiquark, and color-octet state, respectively. The parameter Rs(0) is a fictitious 
"color-octet radial wavefunction at the origin" related to the nonperturbative matrix element 
Hg(A) by H' 8 = (2/37r)|i? 8 (0)| 2 /m^. The square of the amplitude A a for the subprocess 
g* — > cc, summed over final state colors and spins, is 

^AaA} = Q7ra s m 3 c H' 8 (A) (-g a/3 + ^j . (20) 

The q a qp term can be dropped because q a is of order m 2 /q when contracted with the nu- 
merator of the virtual gluon propagator in axial gauge. The expression therefore reduces 
to the polarization sum (—gap + •••) for an on-shell gluon multiplied by 67ra s m 3 H 8 . Di- 
viding by (4ml) 2 fo r the virtual gluon propagator, we obtain the fragmentation probability 
(37ra s /8)Hg/m c . This probability can be identified with the second term in (JJ), integrated 
over z and summed over J = 0, 1,2. This term in the fragmentation function contributes 
only at the endpoint z = 1. We can therefore identify the function ds(z) in ([[]) to be 

d 8 (z) = ^8(l-z). (21) 



The total fragmentation function at leading order in a s is given by the factorization 
formula ([I]), with the color-singlet coefficient given in ( |T4T ) and the color-octet coefficient 
given in (pTj) . To avoid large logarithms of m c /A in the color-singlet coefficient, we can 
choose A = m c . We thus arrive at the final expressions for the fragmentation functions of 
gluon splitting into Xcj to leading order in a s : 

. r KaA2m r ) H'Jm r ) , 
+ 2J+ 1 gi- c> 81 c} 5(1 -z) , 22 
24 m c 

where Qj and Pj(z) given by (p!5|)-(^8|). The choice of the scale [i in the running coupling 
constant is independent of the choice of factorization scale A. We have followed Ref. || in 
choosing /i = 2m c , which is the minimum value of the invariant mass of the virtual gluon. 
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2a 2 s (2m c ) R x 
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m. 



If we wish to use a value for the factorization scale A in ( p2| ) which is significantly smaller 
than m c , we should use the solution ([3]) to the renormalization group equation for H$(A) to 
sum up the leading logarithms of m c /A. 

Rough estimates of the gluon fragmentation contribution to the production of the Xcj 
states at large transverse momentum in any high energy process can be obtained by multi- 
plying the cross sections for producing gluons with transverse momentum larger than 2m c 
by appropriate fragmentation probabilities. Integrating the initial fragmentation functions 
fl22|) over z, we obtain the probabilities 

P g ^ x , pa -Rj ag(2mc)gl + (2 J + 1) TO 'NgsN (23) 
9 XcJ 108m c V ; 24m c ' V ; 

where Rq = 5, Ri = 4, and i?2 = 16. Notice that with the choice A = m c for the factorization 
scale, the color-singlet pieces give rise to negative contributions to the initial fragmentation 
probabilities. Requiring that all the probabilities ( p3|) be positive, we obtain an interesting 
lower bound on H' 8 (m c ): 

H . [mc) > Hl . (24) 

Using E x pa 15 MeV, m c = 1.5 GeV, and a s (2m c ) = 0.26, we find H' s (m c ) > 1.4 MeV. The 
estimate Hg(m c ) w 3 MeV obtained in Ref. |TTJ is consistent with this lower bound. Using 
the value H^(m c ) pa 3 MeV, our estimates for the the initial fragmentation probabilities in 
(p3|) are 0.4- 10~ 4 , 1.8- 10~ 4 , and 2.4- 10~ 4 for Xco, Xa, and Xc2, respectively. The production 
of Xcj states contributes to the inclusive rate for production of the 1 charmonium state 
J/ij) through the radiative decay Xcj — ► J/ift+j- Multiplying the fragmentation probabilities 
given above by the appropriate radiative branching fractions of 0.7%, 27%, and 14%, we find 
that the probability of a J/ip in a gluon jet is approximately 8 x 10~ 5 . This is more than 
an order of magnitude larger than the probability 3 x 10~ 6 for the direct fragmentation of a 
gluon into J /if) that was obtained in Ref. 0. 

The methods used above to calculate the fragmentation functions D g ^ XcJ (z) can also 
be used to calculate the distribution of the transverse momentum p±_ of the Xcj relative to 
the gluon jet. This transverse momentum is related to the the invariant mass s of the gluon 
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jet by s = (p]_ + 4ml)/ z +P±/(1 — z )- For the color-singlet contribution, the s-distribution 
is obtained by integrating over z in (|?|). For the color-octet contribution, the s-distribution 
is a delta-function at s — 4m 2 . Adding these two contributions we obtain 

dP^ = 2a|m^ e(s _ Sm , n(A)) + S (s - 4m») , (25) 

ds 81 s 3 (s-4m;?) 24m c 

ds 27 s 2 (s-4m 2 ) 1 mml JJ + 24m c 1 cj ' 1 j 



^ P ^xc2 4a|r?ygi s 2 + 12smg + 96m 4 iraJI^A) 

^ts~ ~ ~ si S 3 (s _ 4m 2) ^ - wA)) + 5 24mc <H S - 4 ™ c ; > 

(27) 

where s m j n (A) = 4m 2 (1 +A/m c ). Integrating over s, we recover the fragmentation probabil- 
ities given in (|23|) . The cutoff-dependence of the color-singlet terms in (p5|) -(p7|) is cancelled 
by the A-dependence of the parameter H' 8 (A) in the color-octet terms. The color-singlet 
terms in these invariant mass distributions were obtained previously by Hagiwara, Martin 
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and Stirling [TJ]], up to an error of Air in the overall coefficient. They did not include the 
color-octet contributions, so their answers were sensitive to the value of the infrared cutoff 
A. In the region near the lower endpoint s = 4m 2 , the distributions (p5|)-p7|) must of course 
be smeared over an appropriate range in p± before they can be compared with experimental 
data. 

We have calculated the fragmentation functions for gluons to split into P-wave quarko- 
nium states to leading order in a s . The fragmentation functions satisfy a factorization for- 
mula with two nonperturbative parameters H\ and H' 8 which can be determined from other 
processes involving the annihilation and production of P-wave states. These fragmentation 
functions are universal and can be used to calculate the rates for the direct production of P- 
wave states at large transverse momentum in any high energy process. They are also needed 
to calculate the total production rate of the 1 states from the fragmentation mechanism, 
since the P-wave states have significant rates for transitions to the 1 states. The fragmen- 
tation probabilities for g — > Xci and g — > Xc2 were estimated to be on the order of 10~ 4 . This 
is large enough that gluon fragmentation into Xcj should account for a significant fraction 
of the Xcj's that are observed at large transverse momentum in hadron colliders. Fragmen- 



tation into \ c j followed by its radiative decay may also account for a significant fraction of 
the J/ip's that are produced at large pr- 



While this paper was being written, we received a paper by Ma |To| , in which the color- 
singlet term in the fragmentation function for g — > Xci is calculated for longitudinally and 
transversely polarized Xci separately. After summing over polarizations, his result agrees 
with ours except near the endpoint z — 1. In the endpoint region, Ma's fragmentation 
function is sensitive to an infrared cutoff. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grants DE-FG02-91-ER40684 and DE-FG03-91ER40674 and also by 
the Texas National Research Laboratory Commission Grant RGFY93-330. 
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Figure Captions 

1. Feynman diagrams for g* — > cc + g which contribute to the color-singlet term in the 
fragmentation function for g — > Xcj- 

2. Feynman diagram for g* — > cc which contributes to the color-octet term in the frag- 
mentation function for g — > Xcj- 
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